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Three different vortex glass models are studied by examining the energy barrier against vortex 
motion across the system. In the two-dimensional gauge glass this energy barrier is found to increase 
logarithmically with system size which is interpreted as evidence for a low-temperature phase with 
zero resistivity. Associated with the large energy barriers is a breaking of ergodicity which explains 
why the well established results from equilibrium studies could fail. The behavior of the more 
realistic random pinning model is however different with decreasing energy barriers a no finite 
critical temperature. 
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The effect of disorder on the behavior of type-II super- 
conductors in magnetic fields is a subject of enormous 
interest both theoretically and experimentally. Much dis- 
cussion has the last few years focused on the possibility 
of a stable vortex glass phase U and some recent papers 
find evidence from simulations for the existence of such a 
phase in three dimensions [0, |l [|. The necessary ingre- 
dients of a vortex glass model is, however, still an open 
question j|, || and recent results actually suggest that 
two of the popular models, the gauge glass || and the 
random pinning model Q], belong to different universal- 
ity classes. 

In two dimensions, which is the focus of the present 
Letter, the different vortex glass models are commonly 
believed to behave similarly. One of these is the gauge 
glass for which equilibrium |], [?], ||, [| analyses show no 
transition. In the present paper we approach this model 
by examining the low-temperature properties of a vortex 
hopping dynamics. The possibility we explore is that 
ergodicity breaking might lead to different conclusions 
when considering dynamics as compared to equilibrium, 
and we actually find results that strongly suggest a finite- 
temperature transition. Interestingly, the behavior of the 
random pinning model is qualitatively different and this 
points to rich and unexpected behavior of the vortex glass 
models even in two dimensions. 

Results in statistical physics often rely on an assump- 
tion of ergodicity. Only if this assumption is valid it 
is permissible to draw conclusions regarding the physi- 
cally relevant time-averages from the ensemble averages 
fiof . Generally speaking this assumption is expected to 
be valid for models with a smooth energy landscape, 
whereas an energy landscape with valleys separated by 
huge energy barriers, could give rise to a breaking of er- 
godicity. The concept ergodicity breaking presupposes 
some kind of dynamics and in the present Letter we 
examine a vortex hopping dynamics where one vortex 
moves one lattice constant at a time. This is the kind of 
dynamics that is typically used for the study of dynami- 
cal properties in vortex glass models O, [l2|, [l3| . 

The possibility of ergodicity breaking means that the 
standard arguments for the absence of a phase transition 



could fail and this applies both to the rigorous analytical 
argument by Nishimori || and the vanishing of the zero- 
temperature domain wall energy with increasing system 
size [JtI [|, pi. The dynamical studies should, however, 
still be reliable, but they arrive at conflicting conclu- 
sions. Whereas some have reported evidence for a zero 
temperature transition [ pd| others are strongly in favor 
of a transition at a finite temperature flSL [l3| . The dif- 
ferent conclusions seem to be due to different judgments 
regarding the reliability of data at rather low tempera- 
tures. 

In this Letter we present a novel zero-T approach de- 
signed to measure the height of the energy barriers in 
vortex glass models. The quantity in focus is the energy 
barrier against phase slips (defined more precisely below). 
Beside probing the possibility of ergodicity breaking this 
quantity is physically relevant since the resistivity in dy- 
namical simulations is inversely proportional to the time 
between phase slips. A phase slip energy barrier that 
diverges as L — > oo would suggest both breaking of er- 
godicity and a vanishing resistivity (= immobile vortices) 
which in turn implies a superconducting low-temperature 
phase [HJ. 

To determine the phase slip energy barrier we start in 
the ground state and perform an exhaustive search of all 
relevant local moves, in the vortex representation. Con- 
sidering the huge phase space of these models this could 
seem like an impossible task and the only reason why it 
is at all feasible is that it is sufficient to focus on some 
relatively few states with rather low energy. Our main 
result is that the phase slip energy, Vl , in the gauge glass 
model grows logarithmically with system size, a finding 
that we interpret as evidence for the existence of a su- 
perconducting low-temperature phase. 

The Hamiltonian for the random gauge XY model is 

H = Y d uU-e j -A ij -l-T ij -*\, (1) 

(«) v J 

where Bi is the phase at lattice point i on a square lattice 
of size L x L (where L is in units of the lattice constant), 
the sum is over nearest neighbors, and A+j is a randomly 
chosen vector potential in the interval [— r7r,r7r), where 



r is the disorder strength. The gauge glass model cor- 
responds to the case r = 1. In the above Hamiltonian 
we also include the twist variable A = (A x ,A y ) [jl5| to 
allow for phase slips. The unit vector picks or A y 
depending on the direction of the link ij. 

The spin interaction is often taken to be a cosines func- 
tion, but to be able to take advantage of the direct rela- 
tion to the vortex representation we have instead chosen 
U(4>) = J(j> 2 /2. With the frustration given by 

f* = h V X Ar ' 

and fluctuating twists fl5|| the vortex Hamiltonian be- 
comes 

H v = -\ 5> r + / r )G(r - r')(^ + / r <) (2) 

r.r' 

where the q r are unit charges, the factor 1/2 compensates 
for the double counting, and G(r) is the lattice Green's 
function 

G M = f^L] 2 V — ~ 1 

1 L J 2 cos k x + 2 cos k v - 4 ' 

For the discussion of the phase slip energy barrier we 
start by considering this quantity in the pure 2D XY 
model. The ground state is trivially given by q r = 
for all r. To consider vortex transport we create a vor- 
tex pair with a negative vortex at the origin and a pos- 
itive at (1,0), and then move the positive one to (2,0), 
(3,0),. . . (i,0). The twist variable A y will change pro- 
portionally to make the current in the y-direction vanish. 
Since {L,0) in a system with periodic boundary condi- 
tions is the same as the origin we are then back again in 
the ground state but with the twist variable A y = —2tt. 
(We use the term phase slip for this kind of process even 
in the vortex representation.) By these steps we have 
transported unit vorticity across the system. In this sim- 
ple case the energy for a configuration with the positive 
vortex at (x, 0) is given by the lattice Green's function, 
G(x,0), and the phase slip energy barrier for the pure 
2D XY model is the value of this energy at the largest 
separation, Vl = G(L/2,0). 

For the following we note that each elementary move 
may be described by a position r and a unit dipole vector 
P: 

q r — 1 —* q r , put a —1 vortex at r, 
<Zr+p + 1 —* q r +p, put a +1 vortex at r + p. 

To keep track of the transport of vorticity we introduce 
the polarization P. An elementary move always changes 
the polarization, P + p — * P. The vorticity transport 
discussed above starts in the ground state and ends in 
the same configuration but with AP = (L,0). The po- 
larization therefore contains some memory of the steps 
taken and may be used to examine the phase slip. 
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FIG. 1: Illustration of a sequence of elementary moves for 
a realization of a 4 x 4-system, corresponding to transport 
of vorticity across the system, AP = (4,0). The upper left 
corner is the ground state. Charge is symbolized by the shad- 
ing of the squares, black represents —0.67, white represents 
+0.72, and intermediate charges are proportionally grey. The 
digits show the polarization relative to the ground state. The 
charge symbols represent the dipole excited in the next step. 
The bar next to the square shows the energy of each configu- 
ration relative to the ground state. 



For the gauge glass model the sequence of elementary 
moves that corresponds to the lowest possible energy bar- 
rier is in general much more complex than in the pure XY 
model. To illustrate this point we show in Fig. [j] such a 
sequence of moves for a certain disorder realization in a 
4x4 system. In this case six steps are needed to trans- 
port vorticity across the system and return back to the 
ground state. We now describe the method used to find 
such complex paths. 

We first search for the ground state by applying the 
standard spin quench algorithm |lrj| a large number of 
times. The lowest energy state is taken to be the true 
ground state when it has been found A rcp = 10 times. 
No change was found in the results when we instead used 
N Iep = 100 for a set of 10 4 disorder configurations with 
L = 4. We then use an algorithm that is analogous to the 
filling of a (energy) landscape with a liquid rising from 
a source at the lowest position. At each time step it is 
the lowest level accessible by the liquid that is invaded. 
The disorder configurations give rise to different periodic 
energy landscapes in polarization space, and the liquid is 
made to rise until there is a connection between these pe- 
riodic copies. The landscape picture is, however, a great 
over-simplification. For each value of the polarization 
there is a large number of possible configurations with 
different energy. The true picture is therefore more of a 
plumber's nightmare with a large network of pipes criss- 
crossing the polarization space but the intuitive picture 
still holds. With an algorithm that slowly increases the 
level one is guaranteed that the first connection found is 
the lowest possible. 
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In the computer, the algorithm consists of repeating 
the following steps: 

1. Generate AL 2 configurations by applying the 4L 2 
possible dipolc excitations to the current configu- 
ration. 

2. Calculate the energy of each such configuration and 
put them in a sorted list, lowest energy first, to- 
gether with their polarization relative to the ground 
state. 

3. Take the first (lowest energy) configuration from 
this list to be the new current configuration. 

4. If this configuration has already been encountered, 
but with a different polarization such that AP = 
(±£, 0) then we are done. Otherwise, go to step 1. 

To make this algorithm work one also needs a list of 
already used configurations. In step 1 the current con- 
figuration is added to that list and in step 4 each config- 
uration is compared to the list. The main output is the 
phase slip energy barrier which equals the highest energy 
used in the algorithm. The data below are obtained by 
averaging the energy barrier from 3000 to 2 x 10 5 dis- 
order configurations. Much as expected, the number of 
iterations, N, of the above algorithm grows rapidly with 
L, N 0.5expL. In our simulations the spin quench 
algorithm that is used to find the ground state is how- 
ever more time consuming than the determination of the 
energy barrier. 

Figure || shows Vl versus L for several different values 
of r. We first focus on the results for r = 1, the usual 
gauge glass model, shown by solid diamonds in panel 
(a). By plotting the data with log scale on the x axis it 
is found that the points with L > 3 to an excellent ap- 
proximation fall on a straight line with a positive slope. 
We consider this to be strong evidence for an energy bar- 
rier that diverges with increasing lattice size and thereby 
also evidence for the existence of a low-temperature phase 
with immobile vortices. As shown in the figure the scal- 
ing holds down to L = 3. This is the same size as in 
the examinations of the domain wall energy in the gauge 
glass H and it therefore seems likely that the logarith- 
mic increase in Fig. || is the true behavior at large L. 
As further support for this belief we note that the log- 
arithmic size-dependence is not a foreign behavior, but 
is rather built into the model at the outset through the 
logarithmic length-dependence in the vortex interaction. 

For weaker disorder, r < 1, the general feature is the 
same with Vl ~ log L. Considering the change in behav- 
ior as r decreases we find that Vl initially changes very 
slowly. The values of Vl for r = 0.9 and 0.8 (not shown) 
are almost identical to the ones for r = 1. For an anal- 
ysis of the data for r < 0.7 we fit Vl — a(r) + b(r) InL, 
shown by solid lines. When r is reduced below » 0.8, 
a(r) slowly increases but the slope b(r) at first remains 
unchanged, Fig. ||c. The slope actually to a good approx- 
imation remains constant down to r — 0.5 where it starts 
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FIG. 2: Size dependence of the phase slip energy Vl for the 
random gauge XY model. The result for the gauge glass 
[r — 1) are shown by the solid diamonds in panel (a) whereas 
the results for the ordinary 2D XY model (r = 0) are shown 
by crosses in panel (b). We always find Vl ~ log L which sug- 
gests diverging energy barriers against vortex motion and a 
superconducting low-temperature phase. Panel (c) shows the 
slope of the lines in (b). There are error bars corresponding 
to one standard deviation for all the data points, but most of 
them are far too small to be visible. 



to increase and finally approaches b(0) = 2n of the pure 
2D XY model. We believe that the behavior of the slope 
b shown in panel (c) is related to the finding [jlj], [l8| of 
a phase with quasi long range order (without ergodicity 
breaking) for r < r c rj 0.37. Figure Ep is also very sim- 
ilar to the phase diagram for the same model recently 
obtained on the basis of Monte Carlo simulations |Q . 

Since the gauge glass model is an often used model 
of a disordered superconductor in an applied magnetic 
field, the above result would quite surprisingly seem to 
suggest the existence of a dissipation-free low tempera- 
ture phase in disordered thin superconducting films. To 
examine this question in more detail we have also stud- 
ied the phase slip energy in the random pinning model 
introduced in Ref . fTT) , which is meant to be a more real- 
istic model of a disordered superconductor in an applied 
magnetic field. The difference is that the random pinning 
model has a constant non-zero magnetic field at each pla- 
quette whereas the random gauge XY model instead is a 
superconductor with random fields that sum up to zero. 
The Hamiltonian for the random pinning model is 



ql (3) 



The frustration (magnetic field) is here homogeneous and 
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FIG. 3: Size-dependence of Vl for the random pinning model 
with three different values for the frustration. The energy 
barrier decreases with increasing L. 

the disorder is instead included through the pinning po- 
tential v r . We follow Ref. |Tl| and let v r be a random 
variable uniformly distributed between —it and 7t, and re- 
strict the possible values for the vorticity to q r = — 1, 0, 1. 

Figure || shows the size dependence of Vl for the ran- 
dom pinning model with / = 1/3, 1/4, and 1/9. Due to 
the requirement that fL 2 has to be an integer it is only 
possible to get data for a few lattice sizes for each value 
of /. The behavior is here markedly different from what 
we found in the gauge glass; Vl decreases with increasing 
lattice size and we therefore conclude that the phase slip 
energy barrier vanishes in the thermodynamic limit. The 
results are therefore very different from Ref. Jll[ where 
it was concluded from a scaling analysis of the I-V char- 
acteristics that the gauge glass model and the random 
pinning model do belong to the same universality class. 
Interestingly, recent simulations strongly suggest that the 
three-dimensional versions of these two models belong to 
different universality classes || Q . 

We also shortly mention our results for a generalized 
XY spin glass, which is given by the same Hamiltonian as 
the gauge glass but with only two possible values for the 
vector potential, Ay = 0, it. Let s denote the fraction 
of links with Ay = 7r. Due to the chiral symmetry every 
ground state has a chiral mirror image q r >—> —q r that 
is also a ground state. We may therefore define a chiral 
energy barrier V£, as the energy barrier that has to be 
climbed to reach the chirally mirrored ground state. 

The phase slip energy Vl is shown in Fig. 0a. For s > 



Vl is independent of size. For s = 0.05 the data seems to 
go as logL suggesting a critical s c > 0.05, but we cannot 
exclude the possibility that Vl saturates and approaches 
a constant for L > 5. The behavior of the chiral energy 
barrier in Fig. |4(b is different, suggesting V£ ~ log L for all 
s > which means that the system at low temperatures 
and with only local moves is trapped in the part of the 
configuration space with the given chirality. 

To conclude, we have examined the existence of grow- 
ing energy barriers and thereby a possible breaking of 
ergodicity in three different 2D vortex glass models. By 
examining the phase slip energy barrier Vl we found dif- 
ferent behaviors in these three models. In the gauge glass 
Vl increases logarithmically which we take to suggest the 
existence of a low-temperature phase with zero resistivity. 
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FIG. 4: The energy barrier for 2tt phase slip Vl (a) and chiral 
mirroring V£ (b) for the generalized 2tt XY spin glass. For 
the ordinary XY model (s = 0) Vl is shown in Fig. ^| while 
V£ would be identically zero. 



In the more realistic random pinning model Vl instead 
decreases in accordance with a transition at zero tem- 
perature. Yet another behavior is found in the XY spin 
glass where Vl is independent of L. 
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